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Abstract 

We present the full 0{a) electroweak radiative corrections to the double Higgs- 
strahlung process e + e~ — > ZHH . The computation is performed with the help of 
GRACE- loop. After subtraction of the initial state QED radiative corrections, we 
find that the genuine weak corrections in the a-scheme are small for Higgs masses 
and energies where this cross section is largest and is most likely to be studied. 
These corrections decrease with increasing energies attaining about ~ —10% at 
y/s = 1.5TeV. The full 0(a) correction on the other hand is quite large at threshold 
but small at energies around the peak. We also study changes in the shape of the 
invariant mass of the Higgs pair which has been shown to be a good discriminating 
variable for the measurement of the triple Higgs vertex in this reaction. 
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1 Introduction 



One of the primary goals of the first phase of a future e + e~ linear collider [TJ El EH, LC, 
is a detailed precision study of the properties of the Higgs. According to the latest 
indirect precision datajl], one expects a rather light Higgs with a mass below the W pair 
threshold. This is also within the range predicted for the lightest Higgs of the minimal 
super symmetric model(MSSM). The LHC will be able to furnish a few measurements 
on the couplings of the Higgs to fermions and gauge bosons jHj but the most precise 
measurements will be performed in the clean environment of the linear collider. In order 
for these precision measurements to be feasible at the LC, the various observables need 
to be well under control theoretically. Nevertheless, although the branching ratios of 
the Higgs have been computed with great precision it is only during the last year that 
the one-loop electroweak radiative corrections to the main production channel at the 
LC, e + e _ — > uuH, has been achieved^ 13 IB]- Even more recent is the calculation [01 
ITU] ITT] of e + e~ — > tiH which is crucial for the extraction of the important ttH vertex. 
Another important issue concerns the reconstruction of the Higgs potential which is the 
cornerstone of the electroweak symmetry breaking. Although measurements of the Higgs 
mass and its couplings to the fermions and bosons may give some information on this 
potential[12j, a more direct probe is through the measurement of the Higgs self-couplings. 
In particular the tri-linear Higgs self-coupling, HHH, may be accessed in double Higgs 
production. Especially for the light Higgs mass this is a daunting task at the LHC [T3] m]. 
In e + e~ the most promising channel in the first stage of a LC, double Higgs-strahlung, 
e + e~ — > ZHH ^3^2]; seems to be the most promising. The aim of this letter is to give 
the full electroweak correction to this process, thus adding to the precision calculations 
of the Higgs profile at the LC. 

The cross sections for e + e~ — > ZHH are rather small reaching about 0.2fb at centre 
of mass energies around 500GeV and for Higgs masses in the interesting range 110 < 
M H < 130GeV as suggested by supersymmetry. However detailed simulations [T7j in- 
cluding the issue of backgrounds have shown that a precision of about 10% on the total 
cross section can be achieved with a collected luminosity of lab -1 especially if a neural 
networks analysis is performed. Other simulations |18| ITU] have shown that by using some 
discriminating kinematics variables, namely the invariant mass of the HH system, one can 
further improve the extraction of the Higgs self-coupling Xhhh- It is important to point 
out that although the foreseen experimental precision on the cross section is about 10%, 
the knowledge of the complete 0{a) corrections is mandatory for this process. Indeed 
we already know that the radiative corrections to the Higgs self-couplings alone, which 
is just a small subset of the full electroweak corrections, receive a Mf correction and are 
therefore potentially large, with M t being the top mass. The leading correction to the 
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tree-level \ HHH coupling writes* 



A gHg = Ag rgg fl- " JfL ), s 2 w = l-M 2 w /M 2 z . (1) 

Mw and Mz are, respectively, the W mass and Z mass and a the fine structure constant. 
For Mn — 120GeV, this constitutes as much as —10% correction to the tri-linear Higgs 
coupling. Although the contribution from diagrams from the triple Higgs vertex to the 
full process are not the dominant ones, this correction alone can amount to about —5% 
to the total cross section [TH]. 

2 Grace-Loop and the calculation of e + e~ —> ZHH 
2.1 Checks on the one-loop result 

Our computation is performed with the help of GRACE-loop[21 j. This is a code for the 
automatic generation and calculation of the full one-loop electroweak radiative corrections 
in the SM.. The code has successfully reproduced the results of a host of one-loop 2 — > 2 
electroweak processes |21j. GRACE-loop also provided the first results on the full one-loop 
radiative corrections to e + e~ — > vvH OH] and more recently on e + e~ — > ttH 0. Both 
these calculations have now been confirmed by an independent calculation |8] ITU] . For all 
electroweak processes we adopt the on-shell renormalisation scheme according to[7|l2*T|l2*2*]. 
For each process some stringent consistency checks are performed. The results are checked 
by performing three kinds of tests at some random points in phase space. For these tests 
to be passed one works in quadruple precision. Details of how these tests are performed 
are given in[71|^]. Here we only describe the main features of these tests. 
i) We first check the ultraviolet finiteness of the results. This test applies to the whole 
set of the virtual one-loop diagrams. In order to conduct this test we regularise any in- 
frared divergence by giving the photon a fictitious mass (for this calculation we set this 
at A = 10 _21 GeV). In the intermediate step of the symbolic calculation dealing with loop 
integrals (in n-dimension) , we extract the regulator constant Cuv = 1/e — 7^ + log47r, 
n = 4 — 2e and treat this as a parameter. The ultraviolet finiteness test is performed 
by varying the dimensional regularisation parameter Cjjv- This parameter could then be 
set to in further computation. Quantitatively for the process at hand, the ultraviolet 
finiteness test gives a result that is stable over 29 digits when one varies the dimensional 
regularisation parameter Cuv- 



*This can, for example, be extracted from the analysis in |20j . We have also made an independent 
derivation. 
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ii ) The test on the infrared finiteness is performed by including both the loop and the soft 
bremsstrahlung contributions and checking that there is no dependence on the fictitious 
photon mass A. The soft bremsstrahlung part consists of a soft photon contribution where 
the external photon is required to have an energy k® < k c <C Ef,. is the beam energy. 
This part factorises and can be dealt with analytically. For the QED infrared finiteness 
test we also find results that are stable over 23 digits when varying the fictitious photon 
mass A. 

Hi) A crucial test concerns the gauge parameter independence of the results. Gauge 
parameter independence of the result is performed through a set of five gauge fixing 
parameters. For the latter a generalised non-linear gauge fixing conditional] has been 
chosen, 



C-gf = - ieaA^ - igcwpZ^W^ + tw^v + m + ikxsh' 
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,,. (d.Z + tz^(v + eH) X3 y - —(d.Ay . (2) 

The x represents the Goldstone. We take the 't Hooft-Feynman gauge with ^ w = ^ z = 
£a — 1 so that no "longitudinal" term in the gauge propagators contributes. Not only 
this makes the expressions much simpler and avoids unnecessary large cancellations, but 
it also avoids the need for high tensor structures in the loop integrals. The use of the five 
parameters, a, /3, 5, k, e is not redundant as often these parameters check complementary 
sets of diagrams. Let us also point out that when performing this check we keep the full set 
of diagrams including couplings of the Goldstone and Higgs to the electron for example, 
as will be done for the process under consideration. Only at the stage of integrating over 
the phase space do we switch these negligible contributions off. Here, the gauge parameter 
independence checks give results that are stable over 27 digits (or better) when varying 
any of the non-linear gauge fixing parameters. 



2.2 The five point function 

The gauge invariance check is also a very powerful check as concerns the reduction of the 
various tensor integrals down to the scalar integrals of rank N. For iV < 5 we perform 
this reduction in terms of Feynman parameters as detailed in our reviewj^J- As for the 
scalar iV = 3,4 integrals we use the FF package[23j. The implementation of the scalar 
five-point function is done exactly as in our previous paper on e + e~ — > vuH 7 j. On the 
other hand we have improved the algorithm for the reduction of the tensor integrals for 
the five point function. A general five point function is constructed out of tensor five 
point functions and is written as 
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M 



T ( 5) = / dJH N® = G ^... p i dJH l,l v ---l p ^ M <5, (3) 



(2tt)" DoD^-Di J (2vr)" DqDx---D a 



where 



D = / 2 -M 2 D i = {l + s i f -Ml s i = Y / Pn 1 = 1,2,3, 4. (4) 

i=i 

Mj are the internal masses, £>j the set (of linearly independent) incoming momenta, / the 
loop momentum and G^ v "' p a tensor that involves the external momenta, the metric or 
anti-symmetric tensors. The scalar integral has N(l) = 1. 

Because the set Si forms a basis for vectors in 4-dimensional space, introducing the Gram 
matrix = Si.Sj one has the identities 

iT=E«W*i. = J2^(i. Sj ). (5) 

We then rewrite an occurrence of I 2 in N(l) in terms of D . For a term with 1^ we use 
Eq.(j3J) and re-express l.Sj as a difference of Dj and D (and loop momenta independent 
terms). Using this trick one is able to write 

4 

N(l) = Y. E ^)D a + F, 

where F is independent of the loop momentum and corresponds then to the scalar five- 
point function which we treat as in|2]. The terms proportional to D a correspond to a 
tensor four-point function of rank-M — 1 or M — 2, starting from a rank-M in the original 
T®. We use REDUCE for extracting E(l) and F analytically. This new algorithm can 
provide much shorter matrix elements compared with the previous one, which used the 
identity 

instead of Z M in Eq.(jSJ. 
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2.3 Input parameters 

Our input parameters for the calculation of e + e~ — > ZHH are the following. We will start 
by presenting the results of the electroweak corrections in terms of the fine structure con- 
stant in the Thomson limit with a~ l = 137.0359895 and the Z mass Mz = 91.1876 GeV. 
The on-shell renormalisation program, which we have described in detail elsewhere [JI], 
uses M\y as an input. However, the numerical value of My/ is derived through Ar[21j 
with = 1.16639 x 10~ 5 GeV -2 ^. Thus, My/ changes as a function of Mjj. For the lepton 
masses we take m e = 0.510999 MeV, m M = 105.658389 MeV and m T = 1.7771 GeV. For 
the quark masses, beside the top mass M t = 174 GeV, we take the set M u = Md = 63 
MeV, M s = 94 MeV*, M c = 1.5 GeV and M b = 4.7 GeV. With this we find, for example, 
that M w = 80.3766GeV (Ar = 2.549%) for M H = 120 GeV and M w = 80.3477GeV 
(Ar = 2.697%) for M H = 180 GeV. 

As well known, from the direct experimental search of the Higgs boson at LEP2, the 
lower bound of the SM Higgs boson mass is 114.4 GeV[2E]. On the other hand, indirect 
study of the electroweak precision measurement suggests that the upper bound of the 
SM Higgs mass is about 200 GeVjl]. In this paper, we therefore only consider a relatively 
light SM Higgs boson and take the illustrative values M H = 120GeV, M H = 150GeV 
and Mu = 180GeV. In fact already for Mu > 2My/ the process loses its interest for the 
LC as not only the cross section decreases but also because the dominant decay mode of 
the Higgs into W's precludes a precision study [Hj of the triple Higgs coupling. 

2.4 Overview of the Feynman diagrams 

The full set of the Feynman diagrams within the non-linear gauge fixing condition consists 
of 27 tree-level diagrams and as many as 5417 one-loop diagrams for the electroweak 0(a) 
correction to e + e~ — > ZHH, see Fig. ^for a selection of these diagrams. Neglecting the 
electron-Higgs coupling, the set of diagrams still includes 6 tree-level diagrams and 1597 
one- loop diagrams. We define this latter set as the production set. To obtain the results 
of the total cross sections, we use this production set. 

TThe routine we use to calculate Ar has been slightly modified from the one used in our previous 
paper on e + e~ — > vvH jj] to take into account the new theoretical improvements. It reproduces quite 
nicely the approximate formula in [53]. For the QCD coupling, we choose as(Mz) = 0.118. 

*In 9 M s = 92 MeV should read M s = 94 MeV. 
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Graph 16 Graph 24 Graph 40 




produced by GRACEFIG 



Figure 1: A small selection of different classes of loop diagrams contributing to e + e~ — ► 
ZHH . We keep the same graph numbering as that produced by the system. The first three 
graphs represent corrections to the HHH vertex. Graph 720 and Graph 805 are box corrections 
representing Z*,j* — > HHZ . Graph 3369 is a t-channel induced diagram with W + W~ — > ZH* 
as a subprocess. Graph 3417 is a typical pentagon diagram. Graph 5415 is a two-point function 
correction. We do not show the QED graphs as they amount to dressing the vertex e + e~ Z* . 



6 



3 Results 



3.1 Tree- level 

To set the stage and to help understand the behaviour of the QED corrections, let us first 
give a brief summary on the cross section at tree-level. This is shown in Fig. El Especially 
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Figure 2: Total cross section for e + e 
forM H = 120, 150, 180 GeV. 



ZHH as a function of the centre of mass energy 



for the lightest Higgs mass, Mh = 120GeV, there is a sharp increase right after threshold. 
For this particular Higgs mass the peak cross section is about 0.16fb, at higher energies the 
cross section decreases rather sharply. This peak cross section decreases with increasing 
Higgs masses. Note, for further reference, that for energies below ~ 1.2TeV there is a 
rather strong dependence on the Higgs mass even in the restricted range of Higgs masses 
that we are studying. For the measurement of the HHH coupling it is most useful to run 
at the maximum of the cross section. For a total integrated luminosity of lab -1 the la 
statistical error corresponds to about a 2% precision. Thus the theoretical knowledge of 
the cross section at 0.2% is more than sufficient. Expectedly at high energies, the cross 
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section shows little dependence on the Higgs mass, at y/s = 1.5TeV, the cross section is 
about 0.06fb. 



3.2 QED corrections: Slicing vs subtracting 



In our description of the computation of the one-loop virtual electroweak corrections, the 
soft QED bremsstrahlung contribution introduces the soft-photon cut-off parameter k c . 
The k c dependence drops out when calculating the full O(a) by including the hard photon 



part, dan 



da e + e - 



*ZHH-y 



(k° > k e ). With 5y = 5y W + 5y ED the virtual loop correction 



(including weak and QED corrections), the total (integrated) 0(a) correction writes 



CO (a) 



00 (1 + 80(a)) = J dcr (1 + 5 V + S s (k c )) + J da H (k c ), 
da (l + 5 EW ) + J da (S^ ED + 5 s (k c )) + J da H (k c 



°"o(l+<V) 
(T (l + 5 W + 5qed) 



QED 
7 V+S 



(fee) 



(6) 



The second and third terms define our default slicing method. The factorised soft-photon 
correction das writes as 



with duo the tree-level 2 
defined through 



da s = da ■ 5 s (k c ) = da ■ fhhik^ < k c ) (7) 
3 differential cross section and the radiator function $ll is 



LL 



(2vr) 3 2fc0 



where fc 7 (p ± ) are photon (e*) 4-momenta. As a default in GRACE-loop the hard photon 



contribution is performed by the adaptive Monte- Carlo integration BASES |2xj and the 
exact matrix elements are generated by GRACE. We then check that within the Monte- 
Carlo integration errors there is no k c dependence, usually at the per-mil level or even 
better (see below). For practically all 2 — ► 3 processes this is sufficient. 

For the process at hand we have also introduced a second method for the calculation 
of the hard part. A difficulty stems from the large contribution from the collinear re- 
gions, 0(ml/s), (of the hard photon radiation) which cancels a large part of the negative 



contribution from Sy+g 



QED 



+ 5s- In order to have a more stable result we use a 
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subtraction technique [2.8 j which is a variant of the dipole subtraction introduced in 
for QCD and in jSU] for photon radiation. The idea is to add and subtract to dan{k, 
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a function that captures the leading contribution of G?cr#(fc c ) and which is much easier 
to integrate. In our case we use the function da <g> fn, where da is derived from the 
tree-level dao but with kinematics taking into account the radiative photon emission. ® 
stands for the phase-space integration of the radiative photon convoluted with 3-body 
phase-space of the tree-level cross-section. Therefore we explicitly write 

J da (1 + 6 V + 6 S ) + J da ® f LL (k° > k c ) 
+ J (da H (k° > k c ) - da ® f LL (k° > k c )) 

= (7 (1 + 6 W ) + J (da ^ ED + da ® III) + / (da H (k° > k c ) - dd Q ® faitf > k c 



QED QED 
7 A a B 



In a® ED of Eq.(jOJ), the convolution over the photon now includes both hard and soft 
photon emission. The two-dimensional integration over the radiator fn is performed 
with the help of the Good Lattice Point quasi Monte-Carlo method, see for example [HI] 
for a description of this method. This ensures appropriate cancellations of the singularities 
between the loop and bremsstrahlung corrections at each point of (the tree-level) phase- 
space. All integration over the remaining tree-level differential cross section including both 
the virtual loop corrections and photon emission is done with BASES. For a® ED we also 
use BASES adapted to a 2 — > 4 process. Here the subtraction means that all singularities 
are made smooth at each point in phase space. 

For the process at hand, which at tree-level proceeds through s-channel Z-exchange, 
the virtual QED corrections form a gauge invariant set which is all contained in the e + e~Z* 
vertex. The dominant initial state QED virtual and soft bremsstrahlung corrections are 
given by the universal soft photon factor that leads to a relative correction [7j 



2a ( , „ \ i k r 3 



2 



$Ts = ( ^ ~ 1) ln £ + 1 L * + IT - 1 ) , L e = Ms/m 2 e ) , (10) 



7T 



E h 4 6 



where m e is the electron mass, Ey, the beam energy (s = 4E%) and k c is the cut on the 
soft photon energy. 

Although this approach of extracting the full QED correction is the most simple one, 
we have also calculated the full QED corrections separately and subtracted their con- 
tributions from the full 0(a). In order to perform this subtraction, the QED virtual 
corrections are generated by dressing the tree-level diagrams with one-loop photons (the 
photon self-energy is not included in this class). Moreover one needs to include some 
counterterms. One only has to take into account the purely photonic contribution to 



(9) 



9 





C^tree = J d(T 


A 


B 


U QED 


Jda^(k c ) 


fda H (k c ) 


&QED 


[GeV] 


[fb] 


[fb] 


[fb] 


[%} 


[fb] 


[fb] 


[%} 


400 


0.08644 


-0.01636 


0.00007 


-18.84 


-0.11971 


0.10344 


-18.83 


600 


0.16312 


-0.00743 


0.00155 


-3.61 


-0.24122 


0.23535 


-3.60 


800 


0.13631 


-0.00075 


0.00284 


1.53 


-0.21091 


0.21300 


1.53 


1000 


0.10985 


0.00130 


0.00344 


4.32 


-0.17594 


0.18067 


4.31 


1500 


0.06983 


0.00201 


0.00360 


8.02 


-0.11887 


0.12448 


8.04 



Table 1: Comparison between the results of the subtraction method and the slicing 
method as described in the text for totally integrated quantities. Mn = 120GeV and 
k c = 10~ 3 GeV. Sq ED (Sqed) is the percentage total QED correction in the subtraction 
(slicing) method. 



the wave function renormalisation constants of the electron. Performing this more direct 
computation is another test on the system. 

Tabled shows a comparison of the total QED correction in our default slicing method 
and our numerical implementation of the subtraction method. We see that the agreement 
is excellent (at the level of 10~ 4 accuracy). The table also makes clear the advantage 
of the subtraction method in that no large cancellation between the part containing the 
virtual contribution and the rest occurs, moreover both are much smaller than the two 
parts obtained in the slicing method. There, as advertised, a large cancellation takes 
part. We note also that / day^g(k c ) in Tabled is in excellent agreement with the result 
obtained by multiplying the tree-level cross section with the universal factor of Eq. (fTU|) . 
It is important to note that the full QED correction are rather large and negative around 
threshold, moderate around the peak and increase steadily for high energies. Most of this 
can be explained from the behaviour of the tree-level cross section, see Fig. El as a boost 
(sort of a radiative return) towards lower energies. 

3.3 The genuine weak corrections 

We now turn to the genuine weak corrections which are the most interesting from a physics 
point of view. We have already discussed how in this process these corrections can be 
unambiguously defined. The total 0(a) correction and the genuine weak corrections as a 
function of energy for our three representative Higgs masses are shown in Fig. El 

Specialising first to Mh = 120GeV, the genuine weak corrections expressed in the 
a-scheme are very small at threshold and also for energies corresponding to the peak 
cross section. At threshold this is in sharp contrast to the full 0(a) cross section which 
is overwhelmingly dominated by the pure QED correction. For the range of the most 
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Figure 3: The full 0{a) relative correction (top panel) and the relative electroweak cor- 
rection Sw (lower panel) as a function of energy for Mh = 120, 150, 180GeK In addition, 
the genuine weak correction 5 W in the scheme is presented for Mh = 120 GeV (dotted 
line). The insert shows 5w for e + e~ — > ZH for the same three Higgs masses. 
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M H (GeV) 




Acr ( a) (fb) 




S w (%) 


600 GeV 


120 
150 
180 


1.6311(2) x lO" 1 
1.0704(1) x 10" 1 
5.7491(5) x 10 -2 


-0.719(3) x 10~ 3 
-3.987(3) x 10~ 3 
-4.136(3) x 10~ 3 


-0.44(1) 
-3.72(1) 
-7.19(1) 


3.17(1) 
5.14(1) 
7.60(1) 


800 GeV 


120 
150 
180 


1.3631(1) x lO" 1 
1.1221(1) x 10" 1 
8.7604(8) x 10 1 


+4.911(6) x 10~ 3 
-0.167(8) x 10- 3 
-2.917(6) x 10 3 


3.60(1) 
0.14(1) 
-3.33(1) 


2.07(1) 
2.27(1) 
2.45(1) 


1 TeV 


120 
150 
180 


1.0986(1) x 10- 1 
9.7667(9) x 10" 2 
8.5245(8) x 10~ 2 


+3.462(6) x 10" 3 
-0.386(8) x 10" 4 
-3.413(9) x 10~ 3 


3.15(1) 
-0.40(1) 
-4.00(1) 


-1.17(1) 
-1.78(1) 
-2.54(1) 



Table 2: Total 0(a) and genuine weak corrections. 



interesting centre of mass energies where the cross section is largest, the (relative) weak 
correction is never more than 4%. Notice also the "spoon-like" behaviour of the weak 
correction at the lowest energies. This also occurs in e + e~ — > ZH as displayed in the 
insert of Figure OF and has also been observed in the so-called s-channel of e + e~ — > uvH 
(see Fig. 3a ofjlj). This seems to be due to a competition between the bosonic and 
fermionic contributions at energies close to threshold but a more detailed investigation 
is needed to confirm the origin of this common feature. Passed the peak, the weak 
correction steadily turns negative reaching about —10% at 1.5TeV. Comparing to the case 
with M H = 150GeV and M H = 180GeV, we see that, in fact, past = 800GeV, the 
Higgs mass dependence of the genuine electroweak correction is small. Larger differences 
between the 3 Higgs masses considered here only appear for the smallest energies around 
threshold. Like for M H = 120GeV, for both M H = 150 and M H = 180GeV the weak 
correction around the peak is also within about 5% and thus well contained. Note that 
the decrease of the weak correction as the energy increases is faster with the higher Higgs 
mass, a trend which is similar to that inherited from e + e~ — > ZH (see insert of Figure OJ). 
Having subtracted the genuine weak corrections one could also express the corrections 
in the scheme by further extracting the rather large universal weak corrections that 
affect two-point functions through Ar. This defines the genuine weak corrections in the 
scheme as 5^ = 5w — 3Ar. For e + e _ — > uuH this procedure helps absorb a large part 
of the weak corrections. Another advantage is that much of the (large) dependence due to 
the light fermions masses also drops out. However one sees, Fig. El that especially for high 
energies, this scheme fails to properly encapsulates the bulk of the radiative corrections. 
This is akin to what happens in e + e~ — > ZH [32J and e + e~ — > tiH[9\ where the bosonic 

§This calculation of the electroweak correction to e + e~ — > ZH has also been performed by GRACE. 
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weak corrections become important (and negative) at high energies. 
To summarise it is perhaps worth to stress that if the total cross section even at its peak 
value would not be measured better than 10% it would then be difficult to observe the 
genuine weak corrections. One should also observe that the approximate leading top 
mass correction to the Higgs self-couplings as given in Eq. (0) and which has no energy 
dependence can not reproduce the bulk of the weak correction. 



3.4 Corrections to the Mhh distribution 




0.05 



100 200 300 400 500 600 700 800 900 1000 

(GeV) 



Figure 4: da/dM HH for M H = UOGeV at y/s = ITeV. We show the tree-level (full 
curve), the effect of including only the genuine weak corrections (dotted curve) and the 
effect of including the full 0(oi) (dashed curve). 



The Mhh distribution has been shown [T%| IT?!] to be a good discriminator for isolating 
the HHH vertex, taking advantage of the fact that the two final Higgs mimic the decay 
product of a scalar (the virtual Higgs in the HHH vertex) [TH] . It is therefore important 
to enquire how this distribution gets affected by loop corrections. We have chosen as 
an illustrative case, Mh = 120GeV at a centre of mass energy y/s = ITeV which is the 
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same configuration that was studied recently in the simulation of 19] which included, 
at tree-level, the effect of an anomalous triple Higgs coupling. In there it was shown 
that a deviation in this coupling affects primarily the lower end of the Mhh spectrum, 
whereas the higher end is little affected. This is rather similar to the effect of the full 
0(a) correction as shown in Fig. 0] especially at the lower edge of the spectrum. The full 
0(a) correction, which for the integrated cross section reaches about +3%, shows in fact 
a depletion from the higher M HH to the lower M HH values, again an effect due to the 
radiative return. After extracting the QED correction, the distribution including the weak 
corrections is hardly, considering the foreseen precision, noticeable though. Therefore we 
conclude that an anomalous coupling could still be distinguished, if large enough, in this 
distribution provided a proper inclusion of the initial QED corrections is allowed for in 
the experimental simulation. 

4 Conclusions 

We have performed a full one-loop correction to the process e + e~ — > ZHH. This is 
a process which is interesting and of importance mainly because it gives access to the 
measurement of the triple Higgs coupling. Our calculation shows that especially not far 
from threshold the QED corrections are large so that a proper resummation of the initial 
state radiation needs to be performed. However in this energy range the cross section 
are modest and probably not measurable. At energies where the cross section is largest 
on the other hand, the corrections are modest especially for the lightest Higgs mass of 
Mh = 120GeV. This applies also to the genuine weak corrections at around the peak of the 
cross sections. Indeed in these regions these genuine weak corrections do not exceed ~ 5% 
and are therefore below the expected experimental precision. We have also investigated 
how these corrections were distributed as a function of the discriminating variable Mhh, 
the invariant mass of the Higgs pair, having in mind the use of this variable for the 
extraction of the triple Higgs vertex. We find that the genuine weak corrections, contrary 
to the QED corrections, hardly affect the shape of the distribution at least for energies 
where this distribution is to be exploited. 
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Note added 

While finalising this paper a calculation of the same process app eared [33]. We have run 
our program with the same set of input parameters and compared the total 0(a) results, 
see Table. El We find a very good agreement, within 0.1%, for centre of mass energies 
up to y/s = 800GeV for all three Higgs masses, Mh = 115, 150,200GeV. As the energy- 
increases the agreement worsens somehow, at y/s = 1.5TeV the agreement is within 0.3% 
but at y/s = 2TeV the agreement is no better than 0.8%. 



1 / /""-I T T\ 

y/s (GeV) 


71 AT / f~\ T T\ 

M H (GeV) 


C^tree (lb) 


0"O(a) (lb) 


So(a) [%] 




115 


0.17493(2) 


0.1629(2) 


-6.9(1) 






0.17491(2) 


0.16282(2) 


-6.91(1) 


500 


150 


0.071834 (6) 


0.06357(6) 


— 11.50(7) 






0.071830(5) 


0.063529(9) 


-11.59(9) 




200 


0. 49611 (3) ■ 10 J 


0.3329(2) ■ 10 6 


-32.90(4) 






0.49606(4) • 10~ 3 


0.332(3) • 10" 3 


-33.0(6) 




115 


0.14156(3) 


f\ At K-J A / /n \ 

0.1471 (3) 


+3.9(2) 






0.14155(1) 


0.14705(2) 


+3.89(1) 


800 


150 


0.11363(2) 


0.1135(2) 


-0.1(2) 






O.H362(l) 


0.11353(1) 


-0.08(7) 




200 


0.07246(1) 


0.0705(1) 


-2.7(1) 






0.072454(7) 


0.07044(1) 


-2.78(1) 




115 


0.07119(2) 


0.0704 (3) 


-1.1(4) 






0.07118(1) 


0.07058(2) 


-0.85(3) 


1500 


150 


0.06684(2) 


0.0634(2) 


-5.1(3) 






0.06683(1) 


0.06359(2) 


-4.86(3) 




200 


0.06165(1) 


0.0569(2) 


-7.7(3) 






0.061644(6) 


0.05707(2) 


-7.42(3) 




115 


0.05021(1) 


0.0473(2) 


-5.8(4) 






0.05021(1) 


0.04773(2) 


-4.95(4) 


2000 


150 


0.04812(1) 


0.0435(2) 


-9.6(4) 






0.048119(5) 


0.04387(3) 


-8.83(7) 




200 


0.04630(1) 


0.0408(2) 


-11.9(4) 






0.046300(4) 


0.04115(3) 


-11.13(6) 



Table 3: Comparison of results between J3-3J/ and our code. The results of [33] are in italic. 
The Born cross section ow? the full order O(o) corrected cross section (Tola) an d the full 
0(a) relative correction 8o( a ) for various Higgs boson masses and centre of mass energies 
are compared. 
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